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ABSTRACT 

The winter-time transport of geostrophically computed 
angular momentum at 500 mb for the period 1 December 1955 
to 28 February 1950 is obtained for Sault Ste, Marie, 
Michigan. This transport is subjected to a cospectral 
analysis to compute the contribution of harmonic waves with 
periods 1-50 days to the total angular momentum transport. 
Computations are made for tv;o other stations, one east, and 
one south of Sault Ste. Marie for comparison of transport as 
related to general circulation processes. 

The author is de^eply indebted to Professor F. L. Martin, 
Department of Meteorology, U.S. Naval Postgraduate School, 
for his encouragement and invaluable aid in the preparation 
of this manuscript. 
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1. Introduction 

Widger [lO] in his study of the flow of absolute 

angular momentim in the atmosphere stated: 

Although even a complete knowledge of the angular 
momentum transfer in the atmosphere cannot by itself 
furnish a solution of the problem of the general 
circulation, any proposed scheme for this circulation 
should include a means for securing the angular 
momentum flows which are observed. For this reason, 
it is essential that all possible observational knowl- 
edge concerning' these processes be obtained. 

Many authors in recent years, Widger ,• Starr , and 

Mintz [lO, 7, 5 ] for example, have provided studies in 

connection with the relationship of the absolute angular 

momentum transport to the general circulation, or the global 

balance of angular momentiim. In general these studies 

provide quantitative findings of the spatial distribution of 

the total momentum transport and establish the fact that 

there is a mean poleward transport. 

This paper, in hopes of contributing to the knowledge 

of angular momentum transport, has as its purpose a cospect- 

tral analysis of the angular momentum transport for waves of 

periods 1-30 days. Due to the limited time available for 

the study, it was confined to the 500-mb level at a station 

in the mid-latitudes. The primary station employed was 

Sault Ste. Marie, Michigan. 

Several other cospectra are included for other 

stations for the purpose of comparison. However, the 

spectrum for these comparative stations are for periods 

2-30 days. 
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2. Mathematical Formulations 

Widger [lo] has shown that the rate at which the ab- 
solute angular momentum of the zonal wind is increasing in 
a portion of the atmosphere bounded by the earth's surface, 
the top of the atmosphere and extending from the pole to a 
given latitude may be expressed as 

J + Jprdd- f dS 

where dt/is any element of volume within the polar cap just 
described, d 5 an element of the surface of the volume, c/dT 
the projection of an element of the earth's surface area on 
the meridional plane passing through a given point, p is 
the pressure, r the radius of the earth, and the east- 
ward frictional stress acts only on the air at the earth's 
surface . 

In the first term on the right \J^ Is the meridional 
component /v* of the velocity into the polar cap. Since the 
vertical velocities at the earth' s surface and top of the 
atmosphere may both be considered zero, this quantity as a 
whole represents the rate at which angular momentum is 
brought into the region by air motion across the boundary. 
The second and third terms on the right represent, respect- 
ively, the rate at which angular momentum is changed due to 
the torque exerted by pressure differentials across mountain 
ranges and the westv/ard torque , jrT^, , exerted by frictional 
drag at the earth's surface. IIov>rever, since this study 
will consider only the transport through a unit area of the 
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latitude wall at the 500 ~nit) surface, the terms involving 
pressure and friction torques may he neglected. No attempt 
is made herein to determine the left side of equation (1) 
from the computed transport at ’^00 mb. 

The absolute angular momentiim of a gram of air at 
latitude ^ may be written 



taken as an element of the latitude wall , JS > talcing 
the first term on the, right in equation (1), the transport 
of absolute angular momentum may be v/ritten as the sum of 
tv/o terms. 



equal to the rate of advection of "relative" angular 
momentum, and the iX* momentum defined as 



v/hich is equal to the rate of advection of angular momentum 
due to the earth's rotation. Since /L'momentim does not 
lend itself to a cospectral analysis and will not be com- 
puted in this study, the chief emphasis is placed upon the 
average angular momentum transport through a unit area and 
per unit time given by equation (3). Therefore, after 
applying the geostrophic approximation for U and AT in 
equation (3) and dividing by N, one obtains the average 
angular momentum transport 



/^1 m M r c^s -/i * > r, / 



( 2 ) 



Furthermore , since V„-AT and the element of surface may be 




(5) 




(^) 
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Ap 



(5) 




v/here the hydrostatic equation has also been used in (5)* 

At first glance equation (5) appears to he a summation 



apart, and thus a representation of equation (5). However, 
a slightly different interpretation was placed upon equa- 
tipn (5)» in that ZiX represents a grid distance of 500 
nautical miles centered on the station for v/hich the 
computation is being made, and A 2^ and are the contour- 
height increments defining the mean geostrophic “wind com- 
ponents 



in the neighborhood of the station. In addition, Ay was 
taken equal • to AX in equation (6). Thus the interpretation 
of equation (5) actually employed here was to regard it as 
a time-average, over N equal time intervals, of values of 
relative angular momentum transport centered at the station. 
The vertical increment for convenience was taken as Aj? equal 
to one millibar and on the right side of equation (5) 
was scaled down to one centimeter instead of the equivalent 
length of 300 nautical miles used for the grid. The units 
of angular momentum transport in this study will then be 

Equation (5) shov;s that the geostrophic relative mom- 
entum transport depends upon the cross covariance betv;een 



of N terms around the latitude circle at a distance A X 




and /VI ::: 



^ - JL 4 ,g^ 
9 ' f AX 



( 6 ) 
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Uj andA^^. It therefore represents the average geostrophic 
relative angular momentum transported across a section of 
the latitude v;all with grid length ^/iz /cm and vertical 
thickness Af^'= / m 6 • 

The mean cross covariance S(t) of two stationary, 
Gaussian time-series U(t) and at lagT, nay he ex- 

pressed hy the convolution integral 

^(r)- ~p M(t) /or t) d t (7) 

0 

Frequently in meteorological applications, S(r) is 
"evenized" by formings the quantity 



5rr)-f 5 c-t) (8) 

2 

which is independent of the sign of lag T . Henceforth, it 
will be understood that the mean cross covariance has been 
evenized. 

The cross-spectrum (or simply cospectrum) Ch of the 
tv/o series C((t) and is defined as that part of the mean 
cross covariance associated v;ith a unit v/ave band centered 
at wave number h in both series. In symbols may then 
be defined as 



c 






L I A' 

p o? 






(9) 



v/here P is the averaging period and 9^(A) and are 

the Fourier transforms of u(t) and respectively, de- 

fined as in 



W- 



L I m 
T 




- -X w-t 

e 



d t 



« 2ir^ 
—— 



(10) 
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The circiimflex indicates the complex conjugate. 

If one now employs the definition of cross covari- 
ance 5cr) for an evenized function, as contained in equation 
(7)» the Faltung theorem leads to the Fourier integral 



expression 

5 (r)-L /Ywcos (11) 

SO that -SCt') and Ccw) are spoken of as Fourier-transform 
pairs. Equation (11) may be inverted to solve for Ch in 
terms of Also, the resulting inversion can be spec- 

ialized to the case of a finite, but long, averaging period 
r t giving 



rr 



^ J S(r) cos 



I ry>a A 



( 12 ) 



A proof of equation (12) is given, from first principles, 
in Appendix I . 

V/hen discrete data are employed at time intervals 
of A.'t , using a maximum lag of the maximum period, 
or fundamental, that can be resolved is 



- 2. hn (13) 

In equation (12), since T appears as a summand, T will be 
written as 



r = At 



^ 1 , t ( 1 ^) 



with -k. taking on each of the discrete values shown in 
(14). Since the averaging period of equations (12) is 
equation (12) leads to 
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I 



(15) 



^ y\ 






[ 



rvc-'V’s^]^z ^ 



% 



CaS 



17 h-fe ^ 
J 



/ n ^ n) 

In the case of C<j and , the resulting coefficients are 

divided by tv/o, since only one-half bandwidth is available 
at n~0 and m~m , Combining equations (5)» (7)i and (8), 
the cross covariance for A lags is 



C r- 3r .c. as ..^ 

'-^4, 2f"CAxf(N->0 







A2, . o AZ , 



( 16 ) 



sinoe the averaging period of equation (7) is now P~(N’' 

The factor two in the denominator of equation (16) results 
from evenizing. 

The cospectral estimates offered by equation (15) 
are termed "raw estimates". Blackman and Tucey Cl] have 
suggested several smoothing procedures, the simplest of 
which involves application of a moving-weighted average def-’ 
ined by 



Cih - .'2>5*Cn-i '+'.J Cv, + 



2SC„ 



4-1 



(17) 



In the case of there is no higher value to include in 

the average and the \ireighted average is modified so that 

Cm ^ ^ CjM-i Cm (IS) 
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3. General Procedures 

The initial step required establishing what data 
were to be used in making computations. Recalling that the 
computed periods are determined by the relation , 

where ho is the maximiim ninnber of lags at discrete time 
intervals At , it is clear that a half-day lag interval re- 
quires 50 lags to obtain a spectrum for periods 1 - 30 days. 
On the other hand, a one-day sample and lag interval re- 
quires only 15 lags to obtain a spectr\im for periods 2-30 
days. Hov;ever, Blackman and Tukey^l^ state that the best 
results are obtained when the longest lag is about 5 ~ 10 
per cent of the length of the record. This would require 
that the record, if possible, consist of 300 - 600 observ- 
ations when using the 500-mb charts produced at O 3 OO GMT 
and 1500 GMT daily, or half this number when using only one 
series of charts. In order to obtain this number of 
observations, it v;ould require a minimum of 5 months data 
and the spectrum would be influenced by seasonal regimes. 

In order to reduce the data-processing load and inter-annual 
influences as much as possible, it was decided to use the 
observations from the three v/inter months of December 1955 1 
January 1956, and Pebruaiy 1956. Thus the longest lag was 
only 17 per cent of the length of the record for this 
particular season. However, no claim is made that the 
cospectra obtained are necessarily representative of those 
for other winter seasons. 

Data were taken from two sources, depending on the 
station being considered. Since heights of both the O 3 OO 
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GIIT and 15 OO GMT 500-mb level v;ere required at Sault Ste. 

Marie, Michigan, the primary station studied in this paper ^ 
facsimile charts transmitted by the U. S. Weather Bureau 
v;ere utilized. In the case of the other stations used for 
comparison, the smallest period computed was only two days 
and the data used were the 15 OO GMT charts published in 
Daily Series Synoptic Weather Maps [8]. 

From the above-mentioned charts, the geopotential 
height differences were obtained with a 5*0 degree (lati- 
tude) grid. Assuming that normal synoptic smoothing \iras 
applied in preparation of the charts, a grid of this dimen- 
sion should not filter out significant short waves that 
have not already been removed by synoptic smoothing. 

In five cases, where the O 5 OO GMT facsimile charts 
were missing, geopotential height differences were computed 
from the actual observed wind listed in Northern Hemisphere 
Data Tabulations L9j» The observed winds were reduced to 
geostrophic values by subjectively compensating for the 
degree of cyclonic or anticyclonic curvature associated with 
the flow. 

The smoothed cospectra of geostrophic relative 
angular momentiam transport were obtained by using these 
data in equations ( 15 ) > (16) » ( 17 ), and (18), and are graph- 
ically displayed in figs. 5 through 6. All computations 
except for the final smoothing, utilizing equations ( 17 ), (18), 
were accomplished on the GDC - 1604 digital computer. 
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4. Limitations and Sources of Error 

The computational results of this study have been based 
on discrete data using equal intervals Atr . In one series 
of computations, indicated by the solid line in fig. 3? the 
data interval was reduced to 12 hours in order to determine 
the effects of sampling fluctuations and of atmospheric 
oscillations. In all other cospectral analyses, the data 
interval was 24 hours. Throughout the study the maximum lag 
was 15 days. 

The period associated with the smoothed cospectral 
estimate is Z\ry\£it/r) with h = 2^ * * ’ • • ^ ho . On the other hand 

the central frequency of the frequency band is 



~ 2m At 

and the highest frequency resolvable occurs for h = m ; this 
frequency is called the Nyquist frequency fm . The 
cospectral estimate corresponds to the frequency band 



(h- to (h-f j4)y2> mA t 

It should be noted that it is difficult to resolve 
frequencies contributing to the co spectrum at the longer 
periods. For example when h - / , the frequency band spans 
the range Vao cycle / day to cycle / day. 

In order to assess the reality of sharp peaks or 
troughs, Panofsky's criterion for X/y distribution of the 
power spectrum has been followed. Strictly speaking a much 
more rigorous method of assessing significance would have 
been achieved by testing the cospectral correlation 
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function defined by 




Here C>^ is the cospectral estimate for the harmonic, 

and Ph,^ and P ^n^^are the power spectral estimates for UL 



and A/~ at the same harmonic. Hov;ever, the power spectra 
computations were not available in time for the completion of 
this aspect of the work. Hence, the more qualitative ap- 
proach indicated in fig. 1 was attempted with some misgiv- 
ings. It must be emphasized that v/hen the term "95% 
confidence limit" is used, it is only by analogy with corres- 
ponding phraseology to be found in Panofsky and Brier [^4 , 
p. 145] • The actual confidence may be at a somev/hat lo^^^er 



pealc or trough could not have happened randomly (at the 95% 
level of belief) provided the peak or trough in the cospectral 



level. 



Applying the 




distribution to the cospectrum, a 



estimate lies outside the range defined as 




(19) 



where r is the mean cospectral intensity between the peak 
and trough, X"(^s/v Ts the chi-square value at the .05 level 
for V degrees of freedom V = 2 “ V-LhlZ) . Accord- 
ing to Panofsky and Brier p. 
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PERIOD 



Pig. 1. Illustration of the test of 

significance at the 95% level of belief. 



A second source of error which may considerably bias 
the accuracy of cospectral estimates is a phenomenon called 
"aliasing". This phenomenon removes energy from higher 
frequencies than those sampled, and may convert it to lower 
frequencies. Pig. 2 illustrates how aliasing may occur. 




Pig. 2. Illustration of the effects of sampling 

variations in giving an aliased spectrum. 
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Suppose a periodic function, such as that illustrated 
in fig. 2, exists for both uCt) and nr(t) , and consider the 
period of the oscillation shown to he 2 At , If lags are 
talcen at intervals of 3 At , one is investigating the 
cospectral estimate at a period of GAt . If the sampling 
times corresponded to points A, B, and C, one obtains a 
rather large "aliased" cospectral estimate for the period G At. 
Hov/ever, if the sampling points had been E, F, and G, the 
aliased cospectral contribution at period GAt would be 
zero. Note that the oscillation could be a real atmos- 
pheric oscillation or, meteorological "noise" of an obs- 
ervational character. The primary atmospheric oscillations 
of periods shorter than one day are the semi-diurnal (12-hr) 
wave and a 10.5-hr wave, originally described by Pekeris [ 5 ] • 
Hov;ever, the amplitude of these waves is quite small. 

Blackman and Tukey [i] have shown that the frequencies 
which are the principal aliases of frequency f„ are 

, etc., where ^/zAtls the Nyquist freq- 
uency. Hov^ever, the effect of aliasing is most prono\inced 
at frequencies near the Nyquist frequency. Griffith et al 
[^23 suggest that' one way of eliminating aliasing from the 
cospectrum is to delete the cospectral estimates for the 
last 20 to ^0 per cent of the higher frequencies. Hov/ever, 
an experiment was performed, using data for Sault Ste. Karie , 
to assess the extent of aliasing. 

Fig. 5 shows tliree cospectra for Sault Ste. Marie 
based on; 

(a) the 1500 GMT reports at 2^-hr intervals (curve A), 
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(b) the 0300 GMT reports at 24-hr intervals (curve B) , 
and 

(c) both the 0500 GMT and I5OO GMT reports at 12-hr 
intervals (curve C). 

It may be noted from fig* 3 that all cospectral est- 
imates at periods greater than 2, 73 days are in fairly close 
agreement. However, curves A and B deviate rather markedly 
at their Kyquist frequency oC Vz cycle / day. 

In the case of curve C there should be little or no 
aliasing at periods greater than tv;o days according to 
Griffith's criterion, since 50 per cent of the spectrum 
lies between the periods one and two days. On the other 
hand, at all other stations it vdll be assumed, henceforth, 
that all cospectral estimates at periods 2.73 days or larger 
are relatively unaffected by aliasing. 
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5. Interpretation of Cospectra 

Figs. 3 through. 6 present the six co spectra computed 
for this study. In these figures the ordinate is the co- 
spectrum and the abscissa, plotted for convenience on a 
logarithmic scale, is the period. 

In assessing the reality of the various peaks observed 
in fig. 3» the procedure embodied in equation (19) vrill be 
applied only to curve C, as the other two curves do not 
deviate significantly except at a period of 3*75 days, which 
will be discussed separately. In order to compute a mean 
cospectral estimate f®r the peak at a period of 10 days, it 
was assumed that the cosnectral estimate at 30 days is the 
apex of a trough. Inspection of the 95% and 5% fiducial 
limits established from the mean cospectral values on each 
side of the peak and trough, indicate that the peaks and 
trough associated with periods 10 and 6 days, respectively, 
are significant. The pealcs and troughs have only a 5% 
probability of randomly exceeding these limits; therefore 
both are considered to be real. 

W 

In the case of c\irve A, the peak appears not to be 
significantly different when compared to the fiducial limits 
computed for curve C. However, limits computed from mean 
values of curve A, although not shown here in order to avoid 
confusion with those of curve C, show the value at the peak 
on curve B is considerably greater than the limit to the 
left side and slightly greater than the limit to the right 
side. Although the pealcs at 3*75 days are not markedly sig- 
nificant from the limit on 'the right side, there is 
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justification for stating these peaks seem to be real when 
one considers both the left and right side limits. Using 
the significance test discussed above, the reality of any 
minor peaks in fig. 3 at periods of 3 days or less is 
doubtful. 

The curve in fig. 4 is the average cospectrum for 
Caribou during the period 1 December 1955 bo 28 February 
1956. A one-day interval between data points was used to 
compute this cospectrum. An inspection of the fiducial 
limits shown in fig. 4 indicates that the peak at 6 days is 
significantly different from the limits on both the left and 
right side of the peak, and it is therefore considered to be 
real. As to the small peak at approximately 5*3 days, the 
confidence interval shown to the left is quite representative 
for this peak. It can therefore be concluded that this peak 
(at 3*5 days) is probably not real, as it does not lie beyond 
the confidence range. However, it can not be inferred from 
this that the peak is a result of sampling variations, since 
the test does not with certainty disprove the reality of the 
pealc. 

The curve in fig. 5 is also for Caribou, but xiras com- 
puted from data taken during the winter previous to that 
used in fig. 4. The limits in this case show the trough and 
peak at 7*5 and 5 days, respectively, to be real. The trough 
at 3.75 days is questionable, since it is not significantly 
different from the limit to the right. 

One of the more striking features in each of the 
cospectra is the large transport associated with the longest 
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period and the fact this transport has the same sign as the 
total transport. 

In general it is difficult to associate a particular 
period of the cospectrum with a specific space wave, since 
the period is a function of both the wavelength and speed 
of the space wave. It, thereby, seems possible for a trans- 
port at a specific period to have contributions from both 
long and short space v^aves. This would be the case if the 
ratio of their wavelengths is equal to the ratio of their 
speeds. However, assuming that the Rossby wave equation 
gives a valid relationship between the wave speed and wave- 
length for any given mean zonal wind, it may be expected 
that long space v/aves move eastward slowly and short ones 
more rapidly. In view of this assumption, the transport 
associated with a period of thirty days may be attributed 
to the contributions from long space waves. 

Purthermore , large values of transport by the longest 
period were to be expected, for an examination of the 
monthly mean 500-mb charts of December 1955 through 1956 
showed Sault Ste. Marie and Caribou were on the v;e stern 
side of a mean long wave trough, hence the net transport 
would most likely be southward (negative). However, for the 
second cospectrum at Caribou (fig. 5)> the period December 
195 ^ through February 1955 was selected by ascertaining 
through monthly mean 500-mb charts that Caribou was east of 
the trough and a northward (positive) transport could be 
expected. In the case of Athens, Georgia, the main reason 
for selecting this station to provide a comparative 
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cospectrum is the fact it is almost due south of Sault Ste. 
riarie. One could expect these two stations to have a net 
transport of the same sign, but examination of the monthly 
mean 500-mb charts show the mean long wave trough to have a 
southv/est to northeast tilt of sufficient magnitude to place 
the trough west of Athens. As a result the net transport for 
Athens v/as northward, while that of Sault Ste. Marie was 
southward. The indicated convergence of momentum transports 
should then act in such a v;ay as to increase the mean zonal 
v;esterlies. Note also that the only pronounced peak in the 
cospectrum at Athens occurs at a period of 50 days, and that 
the cospectral values diminish v;ith decreasing period. This 
is probably due to the fact that short period waves do not 
penetrate appreciably into this latitude belt. 

Interpretation of significant transports associated 
with shorter periods requires reference to the earlier 
statement, that short v;aves are likely to move rapidly. In 
viev; of this concept the momentum transport evident near 
periods of four days is attributed to waves of approximate 
space-wave numbers 6 to 10, which Saltzmann and Fleisher 
[6] have called "cyclone v^aves". However, the moment\im 
transport peaks at periods near 10 days are presumably due 
to long waves of space-vrave numbers 1 to 5» 
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6. Conclusions 

It nay be concluded, tentatively, from the cospectra 
presented in figs. 3 through 6 and the preceding discussion 
that : 

(a) Haximm transport of geostrophic relative angular 
noHentun is accomplished by disturbances xirith the longer 
periods, 

(b) The pealcs with periods near 10 days are related 
to long space waves. 

(c) The direction of the net transport is hi^ly de- 
pendent on the position of the mean long wave trough 
relative to the observing station. 

(d) Cyclone v;aves of space-wave numbers 6 to 10 may 
contribute to the maximum momentum transports associated 
x;ith periods near 4 days. 

(e) The maximum transports at periods of 50 days may 
be attributed to the circulation index cycle. 
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Fig. 4 Th,e cospectrum of the mean relative angular momentum 
transport at the 500 -mb level over Caribou, Maine for 
the period 1 December 1955 through 28 February 1956. 



Fig. 




5 The cospectrum of- the mean relative angular momentum 
transport at the 500-mb level over Caribou , T'Taine 
for the period 1 December 195^ through 28 February 1955« 
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Fig. 6 The cospectrum of *the mean relative angular momentum 
transport at the 500-mb level over Athens, Georgia 
for the period 1 December 1955 through 28 February 1956. 
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APPENDIX I 



Derivation of the Cospectral Function 

The cross covariance at lagT, of two functions of 
time, is given by 

S(x) 't U(t) . ir('t #■ “) d ■= 

1 0,.)] [ie„ 1 






1 1 /?„ ^ - => ) t)] d 

> »o h I ' 

4 [2»r4y. ) , C ej ^ J 



2dmx 
T . 

27rmr 






c-Z-^r 



The first cosine term integrates to zero for all cases 
m,n ^ 0. The second cosine term integrates to zero for 
all cases m n. When m = n, integration gives 



2 Trio r 
T 



f r 

CO _ _ 

b-a ' 

The evenized cross covariance equals one-half the sum 
of the cross covariance at T and-T , 

* 

\ - 5(<-r) - L V /^/, By, / , /, ,, n zphX ^ 

i^Cr) 2^ ^ 7 y- C O 5 (, /^ S/i ) c oS y 2^ H . Zi h - /, 2, ‘ ' 

h-0 

Defining the unsmoothed cospectrinn as 



Cos()^ ) 
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the Fourier transform of C(h) is obtained as 



S'w- r Ccn) - ^ - dh 

Multiplying both sides by Cas-^ J v and integrating from 
0 to T gives 

S{r)ca5 (t ■*" dT 

4iY hT 

For all n / 0,CoS — integrates to zero, giving 

C(n)- I S(r) COS —^ - dr 

v/here P is the integration limit and T is the fundamental 
peri od. 
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